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MULTI-USER DETECTION IN CDMA SYSTEMS 

FIELD OF THE INVENTION 
[0001] The present invention relates to multi-user detection in Code 
5 Division Multiple Access (CDMA) systems. 

BACKGROUND AND SUMMARY OF THE INVENTION 
[0002] Code Division Multiple Access (CDMA) is based on spread- 
spectrum technology and is a dominant air interface for the proposed modern 

10 3G and 4G wireless networks. The transmitted CDMA signals propagate 
through noisy multipath fading communication channels before arriving at the 
receiver of the user equipment (UE). In contrast to classical single-user 
detection (SUD) computational techniques, which do hot provide the requisite 
performance for modern high data rate applications, conventional multi-user 

15 detection (MUD) approaches require a lot of a-priori information not available to 
the UE. 

[0003] The conventional detection schemes for CDMA signals only 
exploit second order statistics among user codes. Practically, however, the 
underlying user data symbol sequences are in general mutually (near-) 

20 "independent". This is a key assumption, which enables the application of info- 
theoretic learning approaches such as information maximization and minimum 
mutual information to the realm of CDMA. The use of these computational 
techniques is justified since a wide sense stationary slowly fading multipath 
CDMA environment can be conveniently represented as a linear multi-channel 

25 convolution model. The received CDMA signal can be considered as a sum of 
several non-gaussian random variables generated by the linear convolutive 
transformations of statistically (near-) independent component user variables. 
This linear transformation accounts for the user spreading codes, the cell 
scrambling codes (in case of a cellular architecture), multiple channel paths and 

30 slowly fading channel effects. The present invention estimates a linear 
transformation to counteract, as "optimally" as possible, the effects of the 
channel transformation -resulting in the recovery of the original user signals 
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under the constraint of knowing only the user's signature code (and the 
corresponding cell scrambling code for a two stage implementation). 

[0004] Blind Source Recovery (BSR) is the process of estimating the 
original "independent" user-specific symbol sequences independent of, and 
5 even in the absence of, precise system/channel identification. In typical 
downlink signal processing, where many of the system parameters are 
unknown, including the number and type of codes for co-existing users at any 
instant of time, one can use the blind techniques for better estimation of the 
user-specific signals. Alternately, Blind Multi-User Detection (BMUD) 

10 computational techniques, based on the Natural Gradient Blind Source 
Recovery (BSR) techniques in both feedback and feedforward structures, can 
be used. The "quasi-orthogonality" of the spreading codes and the inherent 
"independence" among the various transmitted user symbol sequences form the 
basis of the proposed BMUD methods. The inventive structures and 

15 computational techniques demonstrate promising results as compared to the 
conventional techniques comprising, e.g., Matched Filter (MF), RAKE and the 
LMMSE methods. The inventive computational techniques can be implemented 
either using the batch or the more computationally efficient instantaneous 
update methods. Although batch implementations exhibit better performance, it 

20 is however accompanied by longer latency and require more involved 
implementation structures not suitable for a UE/MS. The remaining text focuses 
on the instantaneous (or on-line) performance of the BSR computational 
techniques, which exceeds the performance of other approaches. However, 
the invention can be easily described in the context of the batch processing. 

25 [0005] This (on-line) detection technique can be easily extended to 

CDMA implementations, using relatively short scrambling codes, but becomes 
impractical in WCDMA downlink using long scrambling codes. In spite of the 
fact that very low bit error rates (BER) can be achieved with the BSR technique 
and the detection process does not even require the knowledge of user's own 

30 signature code, the recovered signal stream is at the symbol level with no 
explicit user identification. Further, inherent sign and permutation ambiguities 
exist in BSR (scaling is not relevant as the recovered streams are typically 
desired to have a constant amplitude (e.g., BPSK, QPSK etc.). User 
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identification in BMUD is not possible unless some preamble or pilot data is 
transmitted periodically. This periodic requirement stems from the dynamic 
nature of the wireless communication scenario where users may dynamically 
enter or exit the system. The environment structure also varies widely due to 
5 the mobility of the MS/UE and the transient in the dynamic environment. 

[0006] With these practical constraints in mind, new computational 
techniques are proposed by an infusion of info-theoretic learning computational 
techniques such as static Blind Source Recovery (BSR) (or Independent 
Component Analysis, ICA) and Principal Component Analysis (PCA) into the 

10 existing structure of a RAKE receiver. The purpose of this additional info- 
theoretic stage is to counter, as best as possible, the unmodeled multiple 
access interference (MAI) and the additive noise contribution of the channel. 
Further, use of a simple info-theoretic stage does not make the receiver 
structure too complex (in fact, it is simpler than most other proposed adaptive 

15 LMMSE implementations. RAKE-PCA uses up to second order statistics, as 
compared to RAKE-BSR, which utilizes higher order statistics. This results in 
slightly simpler update structure for the RAKE-PCA, but the performance of the 
RAKE-BSR is found to be better than RAKE-PCA. Further, assuming the 
score-function for the ICA update law to be chosen properly, the resulting 

20 equalization matrix in case of RAKE-BSR has relatively smaller element values 
(energy) as compared to the corresponding matrix for RAKE-PCA, which can be 
translated to the need of fewer memory bits for storage of coefficients. Lastly, 
both RAKE-BSR and RAKE-PCA use all the available user information, so that 
there are no issues of user identification in this case. The main advantage of 

25 both the adaptive RAKE-BSR and RAKE-PCA computational techniques is the 
improved BER performance for the UE/MS without the need of any additional 
information than what a standard RAKE receiver already has. The proposed 
computational techniques can be applied directly to both generic direct 
sequence (DS-)CDMA and modern multi-cellular 3G (UMTS) and beyond 

30 CDMA systems. The described processes can also be extended to other forms 
of spread spectrum system. 

[0007] Further areas of applicability of the present invention will 
become apparent from the detailed description provided hereinafter. It should 
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be understood that the detailed description and specific examples, while 
indicating the preferred embodiment of the invention, are intended for purposes 
of illustration only and are not intended to limit the scope of the invention. 

5 BRIEF DESCRIPTION OF THE DRAWINGS 

[0008] The present invention will become more fully understood from 

the detailed description and the accompanying drawings, wherein: 

[0009] Figure 1 is a block diagram illustrating a typical signal 

generation model for a QPSK DS-CDMA system; 
10 [0010] Figure 2 is a block diagram illustrating a feedforward demixing 

structure in accordance with a first embodiment of the present invention; 

[0011] Figure 3 is a block diagram illustrating a feedback demixing 

structure in accordance with a second embodiment according to the present 

invention; and 

15 [0012] Figure 4 is a block diagram illustrating a feedback demixing 

structure in accordance with a third embodiment according to the present 
invention. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
20 [0013] The following description of the preferred embodiments is 

merely exemplary in nature and is in no way intended to limit the invention, its 
application, or uses. 

[0014] The present invention includes three embodiments of demixing 
structures providing new MUD detection systems and methods, and two 
25 additional new types of detectors derived using BSR techniques. The new 
MUD detection systems and methods are discussed as Natural Gradient Blind 
Multi-User Detection (BMUD) computational systems and methods. The two 
new detectors are RAKE-Blind Source Recovery (RAKE-BSR) and RAKE- 
Principal Component Analysis (RAKE-PCA) Detectors. These detection 
30 systems, methods, and detectors are discussed with reference to a convenient 
convolutive signal model representation of DS-CDMA systems discussed with 
reference to Figure 1 . 
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[0015] In a typical downlink synchronous DS-CDMA system employed 
for indoor ATM and certain ad-hoc wireless networks, each user's data 10 is 
spread using an individual signature waveform (or spreading code), then the 
data 10 for all users is combined and transmitted over multipath AWGN channel 
5 12 by the Base Station (BS) 14. Each User Equipment (UE) or Mobile station 
(MS) synchronizes itself with the BS using the broadcast synchronization/pilot 
channels; once synchronized, the BS and UE/MS can communicate on the 
traffic channel (comprised of both data and control streams), assuming the data 
transmission to be QPSK, i.e., comprised of two composite data channels 

10 created by a serial-to-parallel (S/P) stage, which are constellated in quadrature. 
At the UE/MS receiver, the received signal is first passed through a chip- 
matched filter and sampled at chip rate. 

[0016] Considering a total of K active users in an L multipath 
environment and N transmitted symbols during the observation frame T F , the 

15 received signal is given by 

N K L-l t ^ (A\ 

K0=Z Z 2 <fi3fr*(»h(tMt-'n-T t ) + n{t) K } 

«=1 k=l 1=0 

where ^is the energy of the n th symbol for the ft'* user, b k (n)e{±l±i} is the 
n th complex symbol for the k th user, fyand v t are the path's gain co-efficients 
20 and delay, respectively. n(t) is the additive noise and s k (t) is the k th user's 
signature code (or spreading sequence) generated by 

*a = (»)/>(' ~ m X ); <x k (m) e {- 1,1}, 0 < m < G - 1 (2) 

a k (m) is a real spreading sequence (i.e., any of the standard CDMA PN codes, 
such as the Gold, Walsh-Hadamard, Kasami sequence, etc.) for the k th user 
25 containing G chips per symbol, i.e., G = T b /T c ,p(t) is a chipping pulse of 
duration T c ., and where being the symbol period. 

[0017] Under the assumption of time-invariance, the model (1) can be 
more compactly written in a vector-matrix format as 

r=HSb + n (3) 
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where, H is a (NG + L-l)xNG multipath propagation co-efficient matrix 
containing the channel coefficients. S is a NGxNK block diagonal matrix with 
the matrix of spreading codes forming the diagonal elements, b is an NK -c/ 
vector containing the user symbols, while n is the (NG + L-i)-d channel noise 
5 vector with covariance matrix £?. The structure of the above defined matrices 
and vectors is given by 

\h 0 o o I 

|i •. : | 

|* w \ o I 

|0 \ ho | 

Lo - 0A w J 
5= S J ,5 = s 2 ■•• 

10 [0018] The compact linear model (3) is useful in deriving the closed 

form expression for linear detectors such as matched filter (MF), linear minimum 
mean squared error (LMMSE) etc. for recovery of the transmitted symbol train 
for a desired user. However, the primary disadvantage of this model is the 
prohibitive dimensions of the constituent matrices, especially with longer frame 

15 durations and larger G, the matrices become excessively large, making this 
model unsuitable for any real-time implementation at UE/MS. 

[0019] Alternately, the signal model can be represented as a linear 
convolutive model, i.e., during the symbol time, the received chip data is 
constituted of the chips corresponding to the currently transmitted symbol, its 

20 delayed multipath components as well as delayed chips from some previously 
transmitted symbols and the channel added noise and artifacts. In this 
formulation, G chips arriving at the UE/MS during the n ih symbol time are 
computed as the sum of the chips from L multipaths of the n th transmitted 
symbol and the multipath components of the previous J-1 symbols 

25 (Vz-l,...,n- J-l), where 
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_ fmax(r L ) 



+1 



(4) 



L /G _ 

and max(r L ) being the maximum chip delay in L multipaths (rounded up). The 
n' h received symbol data can be expressed as 

r„ W =E K {n\J7Jtj ]T h (t )s k (t - nT- r, ) + n n (f ) + 



1=0 



2 Z (» - iW fi *(»-/) W 2 7 '/ ( ? K (* - (n - 7 )T - r, > nT < f < (» + l)T 



7=1 &=1 



1=0 



(5) 



[0020] Under the assumption that max(r L )<G, the above model can 
be expressed just in terms of the current and the preceding symbols. That is, 
the multipaths with delay greater than symbol period either do not exist or are 
weak enough to be ignored. In this case, the output samples of the chip- 
matched filter can be written as: 



k=l L 



1=0 



1=0 



(6) 



where, z^and z kl are G-d early and late code vectors, i.e., 



^[g-^+iJ ••• ^[gJoj^o 

G-r, 



(7) 



(8) 



and is the discretized delay satisfying the constraint 0<r, <T b . Imposing 

time invariant constraints, the multipath slowly fading environment model (6) 
can be represented in the form 

r„=H 0 £„+HA-i+^ (9) 

where b n and b n _ x are the K-d vectors of current and previous symbol for all 

the K users. H 0 and H x are GxK mixing matrices with the structure 

H o=L H o,o H 01 H 0K J,H 1 =LH 1>0 H u • H 1K J 

such that 
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/=0 
1=0 

and the twosome represent the energy of the current and the previous 
symbol respectively. 

5 [0021] The Natural Gradient Blind Multi-User Detection (BMUD) 

system and method is discussed below with reference to Figures 2-4. Blind 
Multi-user Detection (BMUD) is the process to blindly estimate all the user 
symbol sequences directly from the received composite CDMA signal using the 
Blind Source Recovery (BSR) techniques. BSR framework implies recovery of 

10 original signals from environments that may include transient, convolution and 
even non-linearity. The linear BSR computational techniques have been 
developed for linear convolutive mixing environments by the minimization of 
mutual information (e.g., using the Kullback Lieblar Divergence) using the 
natural gradient subject to the structural constraints of a recovery network. The 

15 natural gradient BMUD network can be adapted either in the feedforward or the 
feedback configuration, in which case the proposed BMUD system and method 
adaptively estimates a set of (filter weight) matrices to counter the linear 
convolutive environment model (9). In each case, the input receives either the 
linear convolutive environment model r n or its whitened version (of the 

20 linear convolutive environment model r n ). Parametric matrices 16A-16C W Q - 

W k are adapted to estimate independent user symbols ;y„at an n th instant 

based on the linear convolutive environment model r n or its whitened version 

r w w . A decision stage 18 sums the component mappings to generate the 

output y n to provide the corresponding user symbol estimates b n also at the n th 
25 instant. 

[0022] The justification for BMUD computational techniques is based 
on the convenient convolutive signal model representation of DS-CDMA 
systems, see (9), and the reasonable assumption that the various transmitted 
user symbol sequences are mutually "independent" as they are generated 
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independently of each other. In this framework, both the transmitted sequence 
and the mixing matrices in the model (9) are unknown to the user. The only 
known entity to the user is the self-identification code. Other available prior 
information is the nature of transmitted data, which in standard CDMA systems 
5 is typically assumed to be quaternary sub-guassian distribution, e.g., QPSK 
data distorted by the multipath channel and additive noise. There exists enough 
information to apply the info-theoretic natural gradient Blind Source Recovery 
(BSR) system and method for BMUD in this case. 

[0023] Further assume that the DS-CDMA channel is not over- 

10 saturated and K<G. The proposed BMUD computational techniques do not 
require any pre-whitening of received data. However, in most modern WCDMA, 
G is chosen to be large enough to maximize spreading gain and so as not to 
limit the number of users- in general K < G. Therefore, it is computationally 
advantageous to pre-process received data for dimension reduction to K which 

15 is the actual number of principal independent symbol sequences (or users) in 
the received data. The process of pre-whitening will also remove the second 
order dependence among the received data samples and some of the additive 
noise. The data pre-whitening can be achieved either online using adaptive 
principal component analysis (PCA) computational techniques or it may be 

20 done algebraically by using a large batch (say N samples) of received data, 
i.e., 

* = l>i r 2 — r^, r N ] 

with the correlation matrix 

1 T 
A c = ——-RR r 

c N-l 

25 Then the whitening is achieved using the filtering matrix 

W = D~ V2 V T 
where 

D: represents the /C-dim matrix of principle eigenvalues of the data correlation 
30 matrix A c . 
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V: represents the KxN matrix of principle eigen vectors of the data correlation 
matrix A c . 

[0024] The whitened version of (9) is given by 
rn = W(H 0 b n + Hfa-l +n n ) = H 0 b n + H x b n ^ (12) 
5 where r n w represents the K - d whitened data received during the n th symbol 
time and H^Hj represent the equivalent KxK convolutive mixing matrices for 
the current and the previous symbol (For compatibility, r n is assumed to be N- 

[0025] BMUD computational techniques blindly adapt a set of matrices 
1 0 to output the independent user symbols estimate y n at the n th instant. This is 
followed by a decision stage to interpret, as best as possible, y n and estimate 
the corresponding user symbol estimates b n also at the n th instant. 

k=Hy n ) (13) 
where #/-(•): represents the (nonlinear) decision stage. 
15 [0026] For the feedforward configuration discussed with reference to 

Figure 2, the BMUD stage output is computed as 

yn=W 0 r^+ZW k r^ k (14) 
£=1 

where, in this case, K refers to the number of feedforward coefficient matrices. 
[0027] The update laws for this feedforward structure can be derived. 
20 For the feedforward parametric matrices W Q and W k , the update laws are 

Aw 0 «(r-^yn)y?K 05) 

and 

AW k - (/ - cp{y n )yV )w k - <p(y n )(r™ k f (1 6) 

where <p{) is an element-wise acting score function, and lis a K-d identity 
25 matrix, and k=1,2,... K. 

[0028] For the initialization of the computational technique, W 0 is 
chosen to be either an identity or a diagonally dominantly matrix, while all other 
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matrices W k are initialized to have either random elements with a very small 

variance or as a matrix of all zeros. Note that no matrix inversion is required for 
the feedforward computational technique. 

[0029] In a second embodiment discussed with reference to Figure 3 
5 and hereinafter referred to as feedback configuration I, the output is estimated 
by 

J K \ 
V k=l J 
The update law for this structure using the natural gradient are given for the 
matrix W Q by 

10 AW 0 - -W 0 (l - <P(y n )yn ) (1 8) 

While for the feedback matrices w k , the update law is 

&W k ocW 0 (p(y n )y^ k ) (19) 

for k=1,2,..., K. The matrices in this case are also initialized in a fashion similar 
to the feedforward case. However, note that at least one matrix inversion is 
15 required in this formulation. 

[0030] In a third embodiment discussed with reference to Figure 4 and 
hereinafter referred to as feedback configuration II, the feedback configuration 
implements the feedback structure without the need for any matrix inversion. 
For this feedback configuration II, the BMUD stage output is computed as 

K 

20 yn=W 0 r™-ZW k y n -k (20) 

The update laws for this feedback structure are given by 

AWb«=(/-^(yn)y?K < 21 ) 

AW*oc(r-^y n )y^ (22) 

again, k=1 ,2,...,K, where K here is the maximum number of filter matrices/taps. 
25 [0031] In case the channel is known or can be estimated, the 

performance of the proposed BMUD computational techniques may be adjusted 
by the diagonalization of the absolute value of the global transfer function. The 
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global transfer function presents the combined effect of the complex mixing and 
demixing transfer functions. For the two symbol convolutive models for the 
case of max(T L )<G, the global transfer function for the natural gradient 
computational techniques in the z-domain are given by: 

5 G = G 0 + G 1 z~ 1 + G 2 z~ 2 (23) 

where, for the feedforward configuration 
G 0 =W 0 H 0 =W 0 WH 0 , 

G x = W 0 H 1 + WxHq = W 0 WH 1 + WjWH 0 and (24) 
G 2 =W 1 H 1 =W 1 WH l 
10 while, for the feedback configuration I 

G 0 =Wo 1 H 0 =Wq 1 WH 09 

G l = Wq\Hi - WO = W^iWHi - Wi)and (25) 
G 2 =0 

and for the feedback configuration II 
15 G 0 =Wo#o =W 0 WH 0 , 

G x - W 0 H 1 +W 1 = WqWHi + Wi and (26) 
G 2 -0 

[0032] The proposed BMUD computational techniques as formulated 
result in recovery of the user symbols directly. The computational techniques 

20 can be conveniently applied to DS-CDMA systems using only user-specific 
spreading sequences. They may also be extended to other CDMA systems 
using relatively short scrambling codes, though the dimension of matrices may 
still become large for DSP implementations in a UE/MS. The WCDMA system 
uses long codes in the downlink, making the application of these computational 

25 techniques impractical because of the requisite dimension of the demixing 
network matrices. 

[0033] The RAKE-Blind Source Recovery (RAKE-BSR) and RAKE- 
Principal Component Analysis (RAKE-PCA) Detectors are now described. 
RAKE-BSR and RAKE-PCA are two new proposed adaptive detectors, which 
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utilize the same knowledge as a RAKE receiver. An info-theoretic adaptive 
weighting matrix of dimension GxG is introduced into the RAKE structure, 
which gives a big performance boost to the RAKE receiver. The performance of 
RAKE-BSR/RAKE-PCA exceeds the performance of LMMSE detectors under 
the conditions of high network congestion, imprecise channel estimation, and 
unmodeled inter-cellular interference etc.. The closed form structure of these 
proposed detectors is given by 

$HWH H r for DS-CDMA Systems 

S&C H WH H r for WCDMA Systems ^ 



h URAKE-ICAIPCA 



i 

where W = diag[AA ---A], and Ais the GxG matrix that is adaptively 

10 estimated either using static BSR (ICA) or PCA computational techniques. 

[0034] It is proposed to adapt the matrix A using the natural 
gradient update laws. However, there exist several other methods for ICA/PCA 
and any other suitable method may be used for these adaptations. This blind 
adaptation of the A matrix has several advantages and improves the 

15 performance of the overall equalization process in several ways. Firstly, it can 
dynamically counter artifacts in the estimated channel co-efficients H . 
Secondly, the channel estimation process (as in RAKE receivers) may be 
limited by the structure (such as the number of fingers) and may estimate only a 
few of the dominant channel parameters. The W stage tends to counteract this 

20 anomaly, as best as possible, and provides better performance than LMMSE in 
such cases. Thirdly, this adaptive stage minimizes the effect of the additive 
channel noise, which may have an intricate underlying unmodeled structure. 
Lastly, the natural gradient ICA/PCA computational techniques inherently 
reduce near-far problems by removing any ill conditioning in the signal space for 

25 all the users in the system. This results in all the mobile users in the system to 
have a BER performance similar to the average BER performance of the 
downlink channel. The matrix A is adaptively estimated using the update laws 

A(k + 1) = A(k) + tj k AA(k) (28) 
where 
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AA(fc) = I U-Pty^y^ > A( ^ for static BSR (or ICA) 
\(I-y(k)y(k) H )A(k) forPCA 



(29) 



and <p(-) is a nonlinear score function which depends on the underlying 
distribution structure of the signals involved. For QPSK signal, a suitable score- 
function is 



[0035] Of these proposed RAKE-BSR/RAKE-PCA structures, RAKE- 
BSR exhibits relatively faster and more stable convergence. However, in 
standard CDMA systems the underlying code structure is chosen to be 
"orthogonal", and thus RAKE-PCA may exhibit lower BER solution if the 
channel impairments are linear. Note that in (27), if the channel estimate H is 
either not available or changes very dynamically, the detector can be estimated 
without using the channel estimate and the structure reduces to Matched Filter 
BSR/PCA, i.e., 



The performance of this structure is better than Matched Filter alone, and 
approaches the performance of a RAKE receiver as the underlying matrix 
A converges. 

[0036] The description of the invention is merely exemplary in nature 
and, thus, variations that do not depart from the gist of the invention are 
intended to be within the scope of the invention. Such variations are not to be 
regarded as a departure from the spirit and scope of the invention. 



<Pi (y0 = v t y t - (Xi (tanh(# Re{ yi }) + tanh(# lm{y f })) 



(30) 




for DS-CDMA Systems 
forWCDMA Systems 



(31) 
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